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Abstract 

Generalising the work of Lenard and Bernstein, we introduce a new, 
fully relativistic model to describe collisional plasmas. Like the Fokkcr- 
Planck operator, this equation represents velocity diffusion and conserves 
particle number. However, unlike the Fokker-Planck operator it is linear 
in the distribution function, and so more amenable to a fluid treatment. 
By taking moments, we derive a new fluid model, and demonstrate the 
damping effects of collisions on Langmuir waves. 

1 Introduction 

In the high energy regime in which relativistic effects dominate - such as lab- 
oratory based laser-plasma acceleration [1] - plasmas are commonly described 
by the collisionless Vlasov equation. This approximation is often justified as 
the timescales governing relativistic processes in an underdense plasma are typ- 
ically much shorter than the average time between collisions. However, recent 
advances in high energy density science have increased the demand for descrip- 
tions of plasma dynamics fully incorporating both relativistic and collisional 
effects [2-4]. 

Recent studies of laser-solid interactions [5] and fast ignition fusion [6, 7] 
sought to circumvent this issue by describing a relativistic electron beam inter- 
acting with a nonrelativistic collisional plasma. However, this docs not disguise 
the need for a fully consistent treatment of relativistic collisional plasmas. 

The standard approach to describing plasma collisions is the Fokker-Planck 
equation [8,9], which supplements the Vlasov equation with a term representing 
diffusion in velocity space. This nonlinear integral operator makes the Fokker- 
Planck equation cumbersome to work with in most cases of interest, and does 
not lend itself to the generation of fluid models. As such, it is often replaced by 
simpler model collision operators [10,11]. 

The approach taken by Lenard and Bernstein [12] was to replace the drift 
and diffusion coefficients in the Fokker-Planck equation with simple, distribu- 
tion independent forms. This gave rise to a collision operator retaining three 
crucial features of the full Fokker-Planck description: it represents velocity dif- 
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fusion; it conserves particle number; and it annihilates a Maxwellian equilibrium 
distribution. For these reasons, it has been widely used in the literature [13-15] 

Given the success of the Lenard-Bernstein equation, together with the cur- 
rent interest in relativistic collisional plasmas, it is somewhat surprising that a 
relativistic generalisation has not previously been considered. In this paper, we 
provide such a generalisation, and explore some of its consequences. 

In section 2 we introduce a relativistic collision operator, making full use of 
the geometry of the unit hypcrboloid bundle, and demonstrate that it reduces 
to that of Lenard and Bernstein in the nonrclativistic limit. In section 3 we use 
this operator to derive a fluid model for a relativistic collisional plasma, and 
show how it modifies the dispersion relation for Langmuir waves. 

Throughout, Latin indices run from to 3, Greek indices from 1 to 3, and 
repeated indices are to be summed over their range. 



2 Linear relativistic collision operator 

The Vlasov equation for the electron distribution function, with additional terms 
representing collisions, may be written in the nonrelativistic limit as 

rg-£(*>f+r g). (i) 

Here, q < and m are respectively the charge and (rest) mass of the electron; E 
is the electric and B the magnetic field. In the full Fokkcr-Planck treatment , the 
drift (.F M ) and diffusion (D^) coefficients involve integrals over the 1-particle 
distribution /. 

Seeking a more tractable equation, Lenard and Bernstein [12] replaced the 
drift and diffusion coefficients in |T]) with the simple distribution-independent 
forms 

-> f3v^ (2) 
D"" -► apy™, (3) 

where fi is an effective collision frequency and a is a measure of the equilibrium 
thermal velocity spread. 

These operators retain three important features of the Fokker-Planck oper- 
ator: they represent diffusion in velocity space; they conserve particle number, 

g+V-(nu) = 0, (4) 

where n is the zeroth and nu the first velocity moment of the distribution /; 
and they possess a Maxwellian equilibrium distribution, 

f M =Acxp(-v 2 /2a). (5) 

The Lenard-Bernstein equation, being much simpler than the full Fokker- 
Planck equation, has been used widely [13-15]. Furthermore, as it is lineaiQ, its 
velocity moments provide a well-defined fluid model. Given the current interest 

1 That is, it is linear in /; however, it forms a nonlinear system when coupled to Maxwell's 
equations. 
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in relativistic phenomena in collisional plasmas [2-7], it is natural to seek a 
relativistic generalisation of (0 [3]) . 

The relativistic form of the LHS of |T]) is commonly used in the literature 
on relativistic collisionless plasmas. It has the form 

Lf = ± a(dl_± Ff , a dl\ (6) 
V dx a m dv^ I ' 



where x° = yjl + v 2 and i M = are the 4-velocity coordinates, v 2 — 8^ u v^v u ', 
and F^a = r)* 1 Fb a arc components of the electromagnetic 2-form. rf h are 
components of the spacctime metric tensor, with signature ( — V H — V). 

Using §3§ in (fTJ), the diffusion contribution to the collision operator is pro- 
portional to the Laplace-Beltrami operator on the flat space of nonrelativis- 
tic velocities. The natural generalisation of this to the relativistic case is the 
Laplace-Beltrami operator A on the unit hyperboloid, 

Af = 3^|4 + + v*vA . (7) 



Q v n y J dv^dv v 

See the appendix for a derivation of this expression. 

It is less obvious how to generalise the drift term @. However, the closest 
relativistic operator which combines with ([7]) to annihilate the Jiittner distribu- 
tion 

fj = A'expi-a-Wl + v 2 ), (8) 
(the relativistic generalisation of the Maxwellian) is 

d 



Df = S[^—{fv»). (9) 

Hence a candidate for the relativistic generalisation of the Lenard-Bernstein 
equation HH3]) is 



dxa ro ^J=^Vl + ^(/^) + aA/j. (10) 
2.1 Relativistic covariance 

Although (JTUJ) is relativistic, it is not manifestly covariant. The drift term is 
not invariant under Lorentz transformations, and the diffusion term, being the 
Laplace-Beltrami operator on the unit hyperboloid (not the unit hyperboloid 
bundle) is also not guaranteed to be form-invariant. Only by recasting (|10j) as 
a fully covariant equation can we be sure it is well-defined. 

To make ^ covariant, the Liouville vector L can be written in terms of 
horizontal and vertical lifts of structures on spacetime M. to the tangent bundle 
TM (see [16] for details): 



(—)" - (±F b a —\ 

\ dx a I V m a dx b J 



(11) 



It is important to recognise that, although (JXTJ) is defined on the whole 
of the tangent bundle, we can nevertheless restrict attention to the physically 
meaningful unit hyperboloid bundle 

H = {(x,x) 6 TM : ip(x,x) = 0}, ip = i^ b {x)x a x b + 1. (12) 
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This follows since the Liouville vector field is tangent to H, Lip = 0. 

That ^ is not covariant is to be expected, since by construction the collision 
operator annihilates the Jiittner distribution which itself is not covariant. 
However, the Jiittner distribution can be written covariantly on the tangent 
bundle 

f = A'cxp(-a- 1 Y a x a ), (13) 

at the expense of introducing the 1-form Y = Y a dx a = dx°. 
Expanding the derivative in @, we can write 

Df = 3/3^1 + v 2 f + 0y/l + ■uV^-. (14) 

Both the scalar \fl + v 2 and the vector \fl + v 2 v f± d/dv fl can be expressed in 
terms of the 1-form Y, and maps from spacetime to the unit hypcrboloid. 

Let T,. x : H x — > TAi be the embedding map of the unit hyperboloid over 
x € Ai into the tangent bundle, and define the map 



Then 



t : A X M -► A°TM 
u a dx a H-* u>Yx a . 



(15) 

(16) 
(17) 



where E* is the pull-back by T, x , and jj is the dual with respect to the metric g 
on H x , that is, g(ui^,X) = w{X) for all vectors X on H x . 

Then, defining f x = E*/, the relativistic Lenard-Bernstein equation (Ti"0|) 
can be written 

KLf = 0(3KiY + (dK^f + aA) U (18) 

Since each term in this equation is defined intrinsically, it follows that (fT8|) is 
fully covariant, and hence ([TO]) is well-defined. 

It remains to interpret the 1-form Y . This represents the frame in which the 
equilibrium Jiittner distribution is 'at rest' (that is, its first velocity moment has 
no spatial component) . Given the role of the background ions in determining this 
frame, it is natural to identify Y with the velocity field of the ions, Y = — V\ on . 
The metric dual of a vector X on M. is defined by X(Z) = r}{X, Z) for all 
vectors Z on Ai, rj being the spacetime metric. 

2.2 Nonrelativistic limit 

Having obtained a linear equation for the 1-particlc distribution incorporating 
collisions, we should confirm that it docs indeed reproduce (JTHSJ) in the non- 
relativistic limit. To do this, we introduce a small parameter e, and define the 
scaled velocities 

= ev^. (19) 
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Taking the typical time scale to be of order unity, (|19[) implies that typical 
spatial scales should scale with e, so introduce rescaled coordinates 

x" = x°, x» = ex^. (20) 

In the nonrclativistic limit, we expect a to correspond to the square of the 
thermal velocity spread, and f3 to the collision frequency, which should be the 
inverse of the typical time scale. Hence 

a = e 2 a, (3 = (3. (21) 

Finally, if the electric and magnetic contributions to the Lorentz force are 
to be comparable, components of the electromagnetic 2-form should scale as 

F» V = F%, F\=eF\. (22) 

Then taking the limit s — > 0, keeping circumflcxed quantities constant, (1101) 
becomes 

lh + ~ i ^ + ^ w - *m 0* + ^ w) ■ (23) 

which is indeed ([T}|3|). 



3 Macroscopic fluid models 

Equation pop forms a closed system when coupled to Maxwell's equations, using 
the electric current components 

f = q [ f x *^=L=- qN ^ (24) 

where the ion number current N? on = —ni Qn r) ab Yb, with ni on the (constant) ion 
number density. However, it is not always convenient to work with differential 
equations on the 7-dimcnsional unit hypcrboloid bundle, and solving for the full 
distribution function is often unnecessary. 

Often it can both be more convenient and facilitate interpretation to work 
with a fluid model on spacetime rather than a kinetic model on the unit hy- 
pcrboloid bundle. Such a fluid model may be constructed by taking velocity 
moments of the kinetic equation. Unlike the full Fokker-Planck equation, the 
relativistic Lenard-Bernstein equation (|10p is naturally suited to such a con- 
struction. 

With the usual definition of the number current and the stress-energy tensor 
as respectively the first and second moments of the particle distribution, 

fi a ^£==, S»» = m [fx« x *-^L=, (25) 



vT+i? J VT 

the integral of the relativistic Lenard-Bernstein equation (|TU|) shows that particle 
number is conserved!! 

d a n a = 0. (26) 



2 This equation is valid in an inertial coordinate system. In an arbitrary coordinate system 
V a ". a = 0, with V a the covariant derivative. Similar comments apply to the other fluid 
equations. 
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That the drift term does not contribute a source to the RHS of (J^SJ) is 
evident, as it is the integral of a divergence, and / is assumed to vanish at 
infinity. The same can be seen of the diffusion term, when it is recognised that 
the Laplace-Beltrami operator can be written 

A/ = -yi=-^-(y/d^ g^^-), (27) 

where g is the metric on TL X (see appendix), and the determinant of this metric 
is 

det - 9 = TT^- (28) 

Multiplying (fTU|) by mx b before integrating yields 

d a S ab = -qF b a n a + (d a S ab ) dr . H + (d a S ab ) dm . (29) 

The diffusion contribution to the divergence of the stress-tensor is readily 
calculated using Green's second identity: 

(d a S ab ) diff = map f x b Af 



= ma(3 I Ax b f 



= 3ma(3n b , (30) 

where the result Ai a = 3x a has been used (see appendix for details). 

The contribution of the drift term to the divergence of the stress-energy 
tensor can be integrated by parts (again neglecting boundary terms by assuming 
/ and its derivatives fall off sufficiently rapidly for large v): 

d , „ „, d 3 v 



= /?(7 ? ac y"-7 7 bc y a )5 ac , (3i) 

where Y a = -rf b Yb are components of the 1-form Y = dx° representing the 
velocity of the background ions. 

The full equation for the divergence of the stress-energy tensor then is 

d a S ab = -qF b a n a + f3( V ac Y b ~ ?] bc Y a )S ac + 3maf3n b . (32) 

Note that the term involving the equilibrium temperature a, which originates in 
the curvature of Ti. x , is a purely relativistic effect, vanishing in the nonrelativistic 
limit. 

In practice, it is useful to recast the fluid equations in terms of the average 
velocity field 

u a = \T x n a , (33) 
where h is the zeroth moment of the distribution 
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and the pressure tensor 



P ab = m [ f(x a - u a )(x b - i> b ) f V (35) 
J Vl + v 2 

Then ((3^1) becomes 

b d b u a = — l -d b P ab - ^F\u b - [3(8 a b + u a u b + \p a b )Y b + 3af3u a . (36) 
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Note that the average velocity u a is not unit normalised, but rather its norm 
is related to the trace of the pressure tensor: 

i lab u a u b = -(\ + -Lp\). (37) 
ran 

In common with the fluid equations derived from the collisionless Vlasov 
equation. (f^BI [3"2|) (or equivalently. (PBI [3^1 13^|) ) do not form a closed system, 
but must be supplemented by an equation of state relating S ab and n a (or 
relating P ab , h and u a ). 



3.1 Langmuir waves in a cold plasma 

One commonly used equation of state is the vanishing of the pressure tensor, 
representing a cold plasma. Since a is a measure of the equilibrium thermal 
velocity spread, it should be zero for a cold plasma. Indeed, setting P ab = 
and contracting (|36j) with u a leads to the consistency requirement a — 0. As an 
illustration of the effects of collisions, consider Langmuir waves in a cold plasma. 
Assuming a cold plasma, (|36l |3"T|) become 

V u u = — i u F-/3(Y + Y(u)u), r)(u,u) = -l, (38) 
m 

where u = u a d/dx a , and V is the Levi-Civita connection on M.. The metric 
dual of a vector X is defined by X(Z) = n(X, Z) for all vectors Z on M. ((3"5|) 
then couples to Maxwell's equations as 

d*F = -q-k(nu + n ion Y), dF = 0, (39) 

with n = \/ —rj ab n a n b the electron proper number density. The contribution of 
the collisions to (|3"8|) is proportional the 3-vclocity of the ions as measured in 
the rest frame of the electrons, and may be interpreted as a friction force. 

The simplest solution to (|38 [l3"9"|) . representing electrons at rest and vanishing 
electromagnetic field, is 

n = n ion , u = -Y, F = 0. (40) 

The metric dual of a 1-form A on M. is inverse to that of vectors: r](A, Z) = A(Z) 
for all vectors Z on Ai . 

To see how collisions influence the plasma dynamics, perturb about (|40p . 

u = -Y + eT, n = m on + ev, F = e$, (41) 

and to leading order in e (|38l |3T)|) become 
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V ? T = ^-z p $ - (3T, F(T) = 0, 
m 

d*$ = q*(vY-m on f), d$ = 0. (42) 

Langmuir waves are solutions to (|4"2"|) representing longitudinal oscillations, 
i.e. $ = EdtAdz, T = Xd/dz, with A and ^ functions of (t, z) only. Assuming 
a temporal dependence exp(— iut), (14211 yields the dispersion relation 

« = -<f±^*-(f) a , (43) 

where o; p = y q 2 n- lon /m is the plasma frequency. Thus, the effect of the collisions 
is to downshift the frequency of the oscillations and damp their amplitude. The 
dispersion relation (|43p has previously been derived for high velocity waves 
in [17], using the nonrelativistic Krook collision operator. 

4 Conclusion 

It is becoming increasingly necessary to describe the dynamics of a plasma in 
which both collisional and relativistic effects are important. The full relativistic 
Fokker- Planck equation is for many purposes too cumbersome, so simpler models 
are sought. 

We have presented one such model that retains many of the defining fea- 
tures of the Fokker-Planck equation, but which remains simple enough to yield 
useful information about collective motion. This model may be regarded as a 
relativistic generalisation of the Lenard-Bcrnstcin equation. 

By taking velocity moments of this equation, we have generated a set of 
fluid equations, which when supplemented with equations of state provide a 
description of a relativistic collisional plasma in terms of fields on spacetime. 
As an illustration of this equation, we have shown how collisions modify the 
dispersion relation for Langmuir waves in a cold plasma. 
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A Geometry of the unit hyperboloid 

The spacetime metric r\ = rj^dx 0, (£> dx b can be lifted onto the tangent bundle 
TM. to give the Sasaki metric [16] 

V D = Vab (dx aH ® dx bH + dx aV ® dx bv ^j , (44) 

where if and V are respectively the horizontal and vertical lifts [18]. Using the 
embedding map of the unit hyperboloid over x £ M. into TAA, this gives 
rise to a metric on 7i x : 

g = KV° = (<V - T^)dv^ ® dv", (45) 
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where w M = 5 fiu v l/ . 

The Laplace-Beltrami operator on the unit hyperboloid TL X acting on scalars 
can be written as 

A/ = #d##, (46) 
where # is the Hodge map of g. Expanding this equation, we have 

A/ = §i* d * dv " + d^L* {dv " a * dvV) - (47) 

We can evaluate the terms involving the Hodge map using the orthonormal 
basis and volume 3-form 

e 1 = . dp = , e 2 = pd9, e 3 = p sin 6 dip, #1 = e 1 A e 2 A e 3 , (48) 
\/l + P 2 

and the relations 

it = psin0cos</j, v 2 = psmOsmip, v 3 = pcos9. (49) 

All the terms in (|47|) involving the Hodge map can be calculated explicitly. 
However, due to the isotropy of H x , it is sufficient to find #d#e?u 3 , #(<iv 3 A 
#dv 3 ) and #(dw 1 A#dv 3 ). 

We have: 

cfo 3 = cos 0dp — p sin 9d9 



Vl + p 2 cos 0e* - sin Oe 2 . (50) 



#cfo 3 = v 7 ! + P 2 cos0e 2 A e 3 — sin0e 3 A e 1 



p 2 \f\ + p 2 sin 9 cos 0d0 A dip P sin 2 0d<p A dp. (51) 



d#du 3 = 3p cos — ^=^= sin 0dp A cZ0 A 

= 3u 3 #l. (52) 
From ((50J) and d l A #e' y = S^ftl, it follows that 

#(dw 3 A #dw 3 ) = 1 + p 2 cos 2 = 1 + (u 3 ) 2 , (53) 

and using 

dv 1 = sin cos <^3<ip + p cos cos yci0 — p sin sin (pdip 



— yl + p 2 sin cos ^e 1 + cos cos y;e 2 — sin tpe 3 , (54) 

we have 

#(dw 1 A #ofr; 3 ) = p 2 sin cos cos ip = v 1 v 3 . (55) 

These results can be readily extended, by direct calculation or symmetry argu- 
ments, to obtain 

= 3w", (56) 
#{dv» A #<fo") = 5^ v + u*V. (57) 

Substituting these into (gTJ) yields (J7J). 
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